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Recenty Abrams and Lloyd |lj have proposed a fast algorithm that is based on a nolinear 
evolution of a state of a quantum computer. They have explicitly used the fact that nonlinear 
evolutions in Hilbert spaces do not conserve scalar products of states, and applied a description 
of separated systems taken from Weinberg's nonlinear quantum mechanics. On the other hand it 
is known that violation of orthogonality combined with the Weinberg-type description generates 
unphysical, arbitrarily fast influences between noninteracting systems. It was not therefore clear 
whether the algorithm is fast because arbitrarily fast unphysical effects are involved. In these notes 
I show that this is not the case. I analyze both algorithms proposed by Abrams and Lloyd on 
concrete, simple models of nonlinear evolution. The description I choose is known to be free of the 
unphysical influences (therefore it is not the Weinberg one). I show, in particular, that the correct 
local formalism allows even to simplify the algorithm. 



I. FIRST ALGORITHM 



Step 1. We begin with the state 

= |0i,...,0 n ) 



where the first n qubits correspond to the input and the last qubit represents the output. 
Consider the unitary transformation acting as follows 

C/|0) = i=(|0> + |l)) 

c/|i) = -L(-|o) + |i>) 

Step 2. 

|V[1]) = U® ■ ■■«>£/ ®11V>[0]) 

n 

4e i*i *»>]°> 

* ii...i„=0 

The input constists now of a uniform superposition of all the numbers < n < 2" — 1. 
Step 3. 

\m)=F\m) 

x 1 

= "T^f E \ii,---,in)\f(ii,---,i n )) 

VZ n...i n =0 

where F is some unitary transformation (oracle) that transforms the input into an output; f(ii, . ■ ■ ,i n ) equals 
Step 4- 

\ip[3}) = U' 1 ® . ® U-\ ®l\il>[2]) 

n 

= ^ E + C-1)^ +1 |1X>) ® • • • ® (|0n) + (-ir» +1 |U) ® \f(k, ■ ■ .,in)) 

ii ...iji—O 
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^ E E (-i) (<1+1)Jl+ -" +(< - +1) *'lii,.--.in>|/(i 1 ,...,i»)) (io) 

ji...jn=0 ii...i n =0 
1 1 

- |0i,.--,0 n )|/(*i *„)> 

+ E E (-l) (4l + lbl + - + (l " + 1)j "|jl,...,Jn)l/(n,..^n)) (ID 

{ji-..i»}/{0 I ...0„}»i...»»=Q 

2" - s , 



|o 1) ...,o„)(— 10> + -|1>) 



^ E E (-i) (<1+1)il+ - +(i ™ +1)i "lii.--- ! in)l/(ii> •••,*«)>■ (12) 

{2l-Jn}^{0l...0 n }il...tn=0 



The probability of finding the input in the state 1 1 , . . . , n ) is 

(13 , 

P(s) is a parabola satisfying P(0) = P(2 ra ) = 1 which shows that it has a minimum in s = 2™ _1 . The minimal 
probability of finding the input in the state |0i,...,0„) is therefore P(2™~ 1 ) = 1/2 and it occurs if s is exactly 
one-half of 2 n . 

Probability of finding f{i\, .. ■,i n ) = 1 is s/2™. This intuitively natural result becomes less natural if one tries to 
prove it by using ([l2|). So let's do it explicitly. Let us begin with rewriting Jig ) in the following form 

1 

2, Tl 

ji—jn {ii...i n ;/(ii,...,i n )=l} 



\m) = ^ E E (-i) (ii+lbi+ - +(t " +lb "iji,...,^>ii> 

jl— jr. {»l...i»i/(*l,-..,in)=l} 
4 E E (-l) (ll + 1Wl+ - + (l " + lb "|jl,...,Jn)|0) (W) 
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ji— jr. {ii...i»;/Cii,— ,»»)=0} 

The probability of finding the flag qubit in |1) is 



ji— jr. {ii...i„;/(ii,...,i„)=l} 



(il+l)jl + ...+(ir.+l)jr, 2 



(15) 



= i£| ^ ( _ 1) Hj 1+ ... + Wr.| (16) 

ix—3n {h— i n ;f(h,— ,*n)=l} 

= ^E|(- 1 ) ? ' J +--- + (- 1 ) ?s ' J T ( 1? ) 

where the vectors i r , r — 1 . . . s, are all different (which is essential for the proof) and j = (ji, . . . ,j n ). 

w) = ^E( s +E(- 1 ) (?+?) ' J ) = ^(«-+EB- 1 ) m3 ) - £ ( 18 ) 

because the sum over j vanishes. 

Step 5. We want to distinguish between the cases s = and s > for small s. To do so we are going to use a nonlinear 
dynamics that does not change the "North Pole" 10) but any superposition of |0) with |1) drags to the "South". The 
violation of orthogonality is called, after Mielnik ||], the mobility phenomenon. 

II. MOBILITY FREQUENCY 

Let us first concentrate on a single-qubit system. The first natural guess is something like 



2 



where 



A = ^(lOXOl - |1)(1|) + V / W(|0)(1| + |1)(0| 
and rj is small but nonzero. The solution of (|l9|) for normalized ipo is 

\ilt) = exp [ - ie((^ Wo) ~ (0|A|0)).4i] |Vo) 

e((^o|^|^o) - (0|A|0))tl |^o) - lAsin 



1 cos 



((^o|A^o)-(0|A|0))i]|Vo) 



Assume 



Then 



I'M 



2" - s 



v/(2"-s) 2 + s 2 v /(2«- s )2 + s 2 



II) 



<^o|A|Vo> 



(2 n -s)(0|+s(l| A (2™-s)|0)+s|l) 



(2™-s) 2 + s 2 
(2" - s) 2 r? - s 2 r? + 2(2" - s)s v /T _ ^2 



(2" - s) 2 + s 2 



The mobility frequency is therefore 



(2™ - s) 2 ri - s 2 rj + 2(2™ - sWl - ?? 2 - 77(2" - s) 2 - rys 2 ~2s 2 r] + 2(2™ - s)s-\/l - if 
u e — e — 77; ~ = e — 



(2™ - s) 2 + s 2 
which for 2™ 3> s gives approximately 



(2™ - s) 2 + ,s 2 



syl — ?7 2 es 



On-l 



which makes the algorithm exponentially slow. 
Let us try therefore another nonlinearity: 

i\wi — etanh 



(V#> 



We find 



u> e = e tanh 
= e tanh 
= e tanh 
w e tanh 



<o|A|o) 

(2" - s) 2 + s 2 



(2™ - s) 2 »7 - s 2 r; + 2(2™ - s)sy/l - rf V 
[r](2 n - s) 2 + rjs 2 - (2™ - s) 2 T] + s 2 rj - 2(2™ - s)sy/l - rf 



(2™ - sfrf - s 2 7? 2 + (2™ - s)si]y/l - rf 
2r]s 2 - 2(2" - s) Sy /l - if 



(2™ - s) V - s 2 rf + (2™ - s)sT)y/l - i] 2 
— e tanh 



'V 1 _ V 2 



2 n - 1 rf 



2 n ~ x rf 



(19) 

(20) 

(21) 
(22) 

(23) 

(24) 
(25) 

(26) 

(27) 



(28) 

(29) 
(30) 
(31) 
(32) 



For 77 of the order of 2 (" x ^ 2 one can obtain a reasonable mobility frequency but this requires an exponentially 
precise control over (0|A|0). 
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III. EVOLUTION OF THE ENTIRE SYSTEM 



The discussion given above applies to a single-qubit (flag) subsystem. The entire system that is involved consists 
of n + 1 systems and therefore we arrive at the delicate problem of extending a one-particle nonlinear dynamics to 
more particles. 

The description chosen by Abrams and Lloyd uses the Weinberg prescription. Several comments arc in place here. 
First, it is known that the Weinberg formulation implies a "faster-than-light telegraph". The version of the telegraph 
especially relevant in this context is the one that is based on the mobility effect It is therefore not clear a priori 
to what extent the fact that the algorithm is fast depends on the presence of faster than light effects. Second, the 
Weinberg prescription is meant to describe systems that do not interact. We have two options now. Either we indeed 
want to keep the flag qubit noninteracting with the input (during the nonlinear evolution) or we allow a nonlinear 
evolution which involves the entire quantum computer. If we decide on the first option we should use the Polchinski- 
type description which eliminates the unphysical nonlocal influences, but the nonlinear evolution of the flag qubit 
is determined by its reduced density matrix (the Polchinski-type description was recently formulated for a class of 
equations more general than those considered by Weinberg in ||; its aplication to interacting systems can be found 
for example in jfj). This is the reduced density matrix obtained by the reduction over all 2" states of the input 
subsystem. Physically this kind of evolution occurs if the nonlinearity is active independently of the state of the n 
input qubits. 

But the very idea of the algorithm is to take advantage of the fact that probability of finding the entire input in the 
ground state esceeds 1/4. It is also assumed that one can turn the nonlinearity on and off. It is legitimate therefore 
to contemplate the situation where the nonlinearity is turned on only provided all the input detectors signal 0. 

At this point one might be tempted to act as follows: Take as an initial condition for our nonlinear evolution 
the product state obtained by projecting the entire entangled state on 1 1 , . . . ,0 n ). The problem with this kind of 
approach is that the "projection postulate" of linear quantum mechanics does not have an immediate extension to 
a nonlinear dynamics. There are many reasons for this but I do not want to discuss it here. At this moment it is 
sufficient to know that it is safer to avoid reasonings based on the projection postulate if nonlinearity is involved. 

I propose an alternative formulation. Assume that indeed the nonlinearity is activated only if the input is in the 
ground state. In principle there is no problem with this because all the different combinations of 0's and l's correspond 
to orthogonal vectors in the 2"-dimensionl Hilbert space of the input and there exists, in principle, an analyzer that 
splits a beam of input particles into 2™ different sub-beams. We can place our hypothetical nonlinear medium in front 
of this output of the analyzer that corresponds to the qubinary zero. 

Let us introduce two projectors: 

= |0 1 ,...,0„)(0 1 ,...,0„|®1 (33) 
P = l (n) ®|0)(0| (34) 

Denote by the state of the entire quantum computer, B = l( n ) ® A, and consider the following nonlinear equation 

i m = etanh ( / ' ' '- - - \ . 1 / ; )P^Bm (35) 
1 ' V(*|P(«)B|*) {^\P( n )PBP\^)J 17 y ' 

Both expressions occuring under tanh are time- independent. In particular, for |\E') = Ej i i ^i 1 ...j n i|ii ) ■ ■ ■ i n )\i) 

(VE|p(">P|tt) _ |^o 1 ...o„o| 2 1 , , 

(*|P(«)PPP|*> |* 0l ...o„o| 2 (0|A|0) ri 1 ' 

The other term is constant since the operators under the averages commute with P^B. The term reads explicitly 

(*|PW|*) E fc l* 0l ...0„ fc | 2 



(*|PWB|*) Ew*S 1 ...o B **o 1 ...o B i<fc|A|0 

We know that 



(37) 



2" — s 

* 0l ...o„o = (38) 

^...Onl = £ (39) 

and therefore the mobility frequency is identical to the one obtained for a single qubit description. The explicit 
evolution of the entire entangled state of the quantum computer is finally 
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|* t )= fl-pW +P^coscj' e t-iP^Bsiacj' e tj\^ ) (40) 

For those of the readers who have played a little bit with faster-than-light telegraphs in nonlinear quantum mechanics 
the basis dependence of the evolution may look somewhat suspicious. There is no problem with this, however, since 
the dependence on p(") reflects our experimental configuration: By changing the projector we change the dynamics 
since we simply put the nonlinear device in a different position with respect to the first analyzer. In the faster-than- 
light problem one gets into trouble if such basis-dependent terms are produced at a distance, and this is typical of 
the Weinberg formulation. 

It may be instructive to discuss what would have happened if we did not assume that the nonlinearity is somehow 
activated in a state dependent way. We therefore assume that the flag system does not interact with the input one. 
For this reason we cannot have any dependence on the basis corresponding to the input particles during the nonlinear 
evolution and we use the Polchinski-type extension of the dynamics which looks as follows 

i\m =etanh( \ ' / r l-B (41) 



(42) 



(#|PBP|# 

The first term under tanh is obviously time- independent. The same with the second one which equals 

(*\pm _ Eh...^ i^i...*„oi 2 _ i 

(Hf\PBP\V) E n ...*„ \%,...i n o\ 2 V V 
as before. The solution for the entangled state of our quantum computer is now 

|* t ) = (l(™+ 1 )cosw £ t-iPsinJ) e t)|*o) (43) 



where u> e has to be determined. 

To do so we first compute the reduced density matrix of the flag subsystem 



T*l...n|*><*| 

= ^ E E(- 1 ) (?+?) ' J |/(*1, • ■ • , Q) ■ ■ • , Ol ( 44 ) 



4* 



= E E i-f fa' ■ • ■ ' ^(*!' •••'*«)! 

3 * 

+ E E^^'Wi, • • • . inWii'i, ■ ■ ■ , 4)1 (45) 
2" — s s 

= ^Ho><o| + -|i><i| 

+ ^ E E^ 1 )^'"!/^!, • ■ ■ , in))(/(<i, ■ • • , i' n )\ (46) 
2" — s s 

= — 10)<0| + -|1)<1| (47) 

because the sums over j vanish. The flag susbsystem is therefore in a fully mixed state. Finally 

2™ 1 



= etanh 
= e tanh 



(2 n - 2s)r\ r\ 
s 



(2 n - 1 -s)riJ \2"- 1 ?7 



tanh (2^) <<: ^ (48) 



so it may pay to act with the nonlinearity on a selected subbeam. 

Returning to the question of exponential precision we should note that the nonlinearity I have chosen leads to 
periodic dynamics and for this reason has a vanishing Lyapunov exponent. One could invent a nonlinear equation for 
a two-dimensional dynamics with a positive exponent (cf. Q) but calculations might be less trivial. 
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IV. SECOND ALGORITHM 



The first three steps are identical to the previous ones 
Step 4- We begin with the result of the third step 

1 



|#2]> = -^= \h,...,in)\f(h,...,in)) 



(49) 



i 1 ...i n —0 

We assume that f(n) — 1 for at most one n (s equals or 1). The state fl49| ) can be written as 

l 



|^[2]} = -^= J2 \0l,i2,...,in)\f(0l,i2,...,in)) 
I 1 



i 2 ...i„=0 

Let us note that with very high probability the state is 

l 



M2]> = -L £ |Oi,i 2 ,...,^>!o> 

v i 2 ...i„=0 

I 1 

+ 7^ E |li^2,..-,i„)|0) 

With much smaller probability it is either 



l 2 ...l„=0 



\m) 



! 1 

;= ^ |0i,i 2 ,---,i«)|l) 

i 2 ...i„=0 

j 1 

7= l 1 !:^, ■ ■ ■ , ^>|0> 



i 2 ...i„=0 



! 1 

,= £ |(M 2 ,...,i n )|0) 

i 2 ...i„=0 

j 1 

— 5^ |ii,j 2 ,..., in )|i) 



i 2 ...i„=0 



and is never in the form 



1 1 

\m) = ~7= E I0l^2,..., in>|l> 
v i 2 ...i„=0 

! 1 

7^ E |ll.*2,..-,*n)|l> 



i 2 ...i n =0 



(50) 



(51) 



(52) 



(53) 



(54) 



since this would mean there are two different numbers satisfying f(n) = 1 which contradicts our assumption. The 
idea of the algorithm is to use a nonlinearity that leaves ([H]) unchanged but (|5^) and ( |53| ) transforms, respectively, 
into 



l« = 4= E I0l,i2,...,i„)|l) 
v j 2 ...i„=0 

! 1 

-= J] i 2 , . . 



(55) 



i 2 ...i„=0 



G 



and 



i« = 4= E |oi,i 2 ,...,i„)|i) 

i 1 

-= 2 |li,*2,...,U|l) (56) 



1 

2™ ^ 

i 2 ...i„=0 

One should be aware of the fact such transformations are in fact impossible within the nonlinear Schrdinger equation 
framework since one cannot merge two different vectors into a single one if the dynamics is reversible and first order in 
time. However one can do this with arbitrary acuracy as can be clearly seen from the preceding examples. The more 
serious problem is that using exactly the same trick it was shown in Q that this kind of evolution leads to influences 
between separated systems (here the flag system would influence the first qubit). In the Weinberg description this 
leads to a contradiction when one obtains this kind of behavior assuming simultaneously that the subsystems do not 
interact. 

We have again two possibilities. We can either assume some form of interaction between the subsystems, or take 
a correct (n + l)-particle extension of a nonlinear dynamics of the flag subsystem assuming no interactions between 
different subsystems. 

Before launching into a more detailed analysis let us first illustrate the Abrams-Lloyd idea on a simple example. 
Take n = 3 and /(HO) = 1. The oracle produces 

^[1000) |0) + 
|001) |0) + 
|010) |0) + 
|011)|0) + 
1 100)|0) + 
|101)|0) + 
|110)|1) + 

|111)|0)] (57) 

The nonlinearity now looks at the second and the third input slots and sees the above kets as the following pairs 

1 



= [|000)|0) + 1 100) |0> 



2V2 

|001)|0) + 1 101) |0> + 
|010)|0) + |110)|1> + 

|011>|0) + |111)|0>] (58) 

Now it scans each of the rows and does not do anything when two flag 0's occur, but when it notices one and one 1 
it changes to 1. So after this step we get 

' [|000)|0) + |100)|0) + 



2V2 

|001)|0) + |101)|0) 
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|010)|1) + |110)|1) + 

|011>|0> + |111>|0>] (59) 

Now the nonlinearity looks at the first and the third slots and sees the kets as the following pairs 

1 



:[|000)|0) + |010)|1) 



2V2 [ 

|001)|0) + |011)|0) + 

|100)|0) + |110)|1) + 

|101)|0) + |111)|0)] (60) 



It again behaves as before and what we get after this step looks as follows 

_L[|ooo)|i) + |oio)|i) + 

|001)|0> + |011)|0) + 
|100)|1) + |110)|1) + 

|101)|0) + |111)|0)] (61) 

Finally our nonlinearity looks at the first and the second slots and the state regroups in the following way 

1 



:[|000)|1) + |001)|0) + 



2V2 1 

|010)|1) + |011)|0) + 

|100)|1) + |101)|0) + 

|110)|1) + |111)|0)] (62) 
Now each row contains one 1 and in the final move all flag 0's are switched to l's and the state partly disentangles: 

-±= [|000) + |001) + 
|010) + |011) + 
|100) + |101) + 

|H0) + |111)]|1) (63) 

Of course, in case s — the entire state does not change during the operation and a measurement on the flag qubit 
gives with certainty. 

One can try to implement such an evolution in terms of a Schrodinger-type dynamics. Let us note that the above 
procedure is somewhat artificial. Once we agree that the nonlinearity can somehow globaly and simultaneously 
recognize the states of all the qubits the optimal strategy would be to choose a nonlinear evolution which changes all 
flag 0's into l's if at least one 1 has been "seen". 
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As has been already said this kind of dynamics is unacceptable if one wants to apply the nonlinear evolution locally 
only to the flag qubit. Let us proceed therefore differently and apply the Polchinski-type description. Begin with the 
nonlinear 1-particle equation 

i\fy=etaxih(a{i)\A-ril\i)))A\i)) (64) 

where 77 and A are the same as before but a is a very large number. For \ip) = |0) the expression under tanh vanishes. 
For a small admixture of |1) and sufficiently large a the mobility with a nonzero frequency begins and an arbitrary 
small amount of |1) can be sufficiently amplified. The extension to the entire quantum computer is 



The solution is 



with 



where 



so 



= etanh (a<tf|l<"> <8> (A-r;l)|tf))l< n > <8> (65) 

|*t) = (l (n+1) cosw £ , Q i - iB sin w e>a t) |* ) (66) 

w £ , Q = etanh (a<*|l< n > ® (A - 

= etanh (aTr p(A - 77I)) (67) 



Tr P (A - ,1) = 2-Tr ( 2 " ~ * ° ) ( ^JL_ ) = -2-2,, (68) 



w £jQ = etanh (^^) (69) 
Now we can explicitly calculate the average of ct 3 = |0)(0| — |1) (1| at the flag subsystem: 

<g> a 3 |* t ) = (* |(l ( " +1) cosw £ , a f + iBsinu; £ia t)l ( ™ ) <g> cr 3 (l (n+1) cosu; £ia t - iBsinw^t) |* ) 
= (*o|cos 2 w £ , Q tl^ ®a 3 |*o) 
+ (* |sin 2 w eia tSl( n )®<7 3 B|*o> 
+i(*o| sinw £ia icosw e , a iBl (n) <g> cr 3 |* ) 
-i(*o| sinw £;Q icosw £iQ il (Tl) <g> a 3 B\^ ) 
2™ — 2s 

= cos 2 uj e . a t + sin 2 w £ . Q t(*o|l (n) O ^cr 3 A|* ) 

z™ 

+ ^sin2w £ , Q t(*o|l (n) ® [A^a]|*o> 
2«— 1 _ s 

= cos 2 uj^ a t +sin 2 w e , a t(* |l (n) ® ^cr 3 A|*o) 



+ _ A /1_ 7? 2 s in2 W£iQ t(* |l(") ® [cti,(T3]|*o> (70) 



Now 



.,2 



?7 2 —77 



1 „ / 2" - * \ / 2t7 2 -1 2r /v /l _ 7f 



2» Tr V 5 y ^ 277v/l- 77 2 l-2r/ 2 



i-(2»-2s)(2^-l) (71) 



2 
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(* t |lW®oa|*t) = 2n _! cos 2 ^+^ ^ sin 2 w e ,at 

2"-i _ s 2™~ 1 — s 

= x cos 2w^ a t + 2rf x sin 2 u} £jCt t (72) 

For s = the average is constant in time and equals 1. For s ^ and rj 2 « and sufficiently large a it oscillates with 
<^e,Q ~ £• This kind of algorithm cannot distinguish between different nonzero values of s, but clearly distinguishes 
between s — and s ^ in a way that is insensitive to small fluctuations of the parameters. 

What is important, such an algorithm is obtained by applying the nolinear evolution only to the flag qubit. This 
is done in a fully local way. We conclude that nonlinear quantum evolutions can lead to fast algorithms and the fact 
that they are fast does not follow from unphysical faster-than-light effects. 
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